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Cerenkov angle and charge reconstruction with the RICH detector of the 
AMS experiment 
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a LIP-Av. Elias Garcia, 14, 1000-Lisboa, Portugal 
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The Alpha Magnetic Spectrometer (AMS) experiment to be installed on the International Space Station (ISS) 
will be equipped with a proximity focusing Ring Imaging Cerenkov (RICH) detector, for measurements of particle 
electric charge and velocity. In this note, two possible methods for reconstructing the Cerenkov angle and the 
electric charge with the RICH, are discussed. A Likelihood method for the Cerenkov angle reconstruction was 
applied leading to a velocity determination for protons with a resolution of around 0.1%. The existence of a large 
fraction of background photons which can vary from event to event, implied a charge reconstruction method based 
on an overall efficiency estimation on an event-by-event basis. 



1. The AMS experiment 

The Alpha Magnetic Spectrometer (AMS) 
is a precision spectrometer to be installed by 2005 
in the International Space Station (ISS), where it 
will operate for a period of three years. Its main 
goals are the search for cosmic anti-matter, the 
search for dark matter and the measurement of 
the relative abundance of elements and isotopes 
in primary cosmic rays. 

The future installation of AMS in the ISS was 
preceded by a 10 days engineering test flight 
aboard the Space Shuttle Discovery in June 1998, 
at a mean altitude of 370 km. Although the pur- 
pose of this experimental flight was the test of the 
spectrometer design principles, about 100 million 
events were collected enabling precise measure- 
ments of the spectra of high energy protons, elec- 
trons, positrons and helium nuclei [Q]. 

The AMS spectrometer capabilities were ex- 
tended with respect to those of the experimen- 
tal flight, through the inclusion of new subdetec- 
tor systems and the completion of others. The 
spectrometer design includes a superconducting 
magnet, a Time-of-Flight system (TOF), a Silicon 
Tracker, Veto Counters, a Transition Radiation 
Detector (TRD), an Electromagnetic Calorimeter 



(ECAL) and a Ring Imaging Cerenkov Detector 
(RICH). It will be capable of measuring the rigi- 
dity (R = pc/\Z\e), the charge (Z), the velocity 
(P) and the energy (E) of cosmic rays within a ge- 
ometrical acceptance of ~ 0.5 m? .sr. The track- 
ing system, with a cylindrical shape, is made of 
8 double sided silicon planes embedded inside a 
magnetic field of about 0.9 Tesla and will pro- 
vide both charge measurements and momentum 
measurements with a resolution Ap/p of at most 
2% up to 100 GeV/c/nucleon. The TOF system, 
made of four scintillator planes placed at the mag- 
net end-caps, will provide a fast trigger, charge 
and velocity measurements for charged particles 
as well as information on their incidence direc- 
tion. On the top of the spectrometer, the TRD 
will discriminate between leptons and hadrons, 
and, on the bottom, the ECAL will contribute to 
the e/p separation and will measure the energy of 
the gamma rays crossing AMS. Figure [l] shows a 
schematic view of the AMS spectrometer. 

The RICH detector will operate between the 
TOF and the ECAL subdetectors. It was de- 
signed to measure the velocity of singly charged 
particles with a resolution A/3//3 of 0.1%, to ex- 
tend the electric charge separation up to the Iron 
element and to contribute to the albedo rejection. 
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Figure 1. A whole view of the AMS Spectrometer. 



Its acceptance is of ~ 0.4 m 2 .sr, that is, around 
80% of the AMS acceptance. The RICH is a pro- 
ximity focusing detector with a low refractive in- 
dex radiator (aerogel) on the top and a pixelized 
photomultiplier matrix on the bottom, where the 
radiated Cerenkov photons are collected. The ac- 
tive pixel size is of 8.5 mm. A conical shaped 
mirror surrounds the whole set, increasing the de- 
tector reconstruction efficiency. Constraints on 
the amount of heterogeneous matter in front of 
the downstream electromagnetic calorimeter have 
imposed a large non-active (empty) readout area 
in the detection plane. For a more detailed des- 
cription of the detector see reference || in these 
proceedings. Figure || shows a view of the RICH 
detector and of its dimensions. 

2. Velocity reconstruction 

A charged particle crossing the RICH radiator 
material of refractive index n, emits photons if its 
velocity (/?) is larger than the velocity of light in 
that medium. The aperture angle of the emitted 
photons with respect to the radiating particle is 
known as the Cerenkov angle, C , and it is given 



by 1: 

cos6> c = -J- (l) 
p n 

It follows that the velocity of the particle, j3, is 
straightforwardly derived from the Cerenkov an- 
gle reconstruction. 

The emitted photons can suffer interactions in 
the aerogel radiator (Rayleigh scattering, absorp- 
tion), can be reflected or absorbed on the mirror 
surface and can fall on the active area composed 
of solid light guides on top of the photomulti- 
pliers. As a consequence, the reconstruction of 
the Cerenkov angle has to deal with two kinds 
of photons; those which are only slightly devi- 
ated from the expected photon pattern due to 
the pixel granularity, radiator thickness and chro- 
maticity effects, and, those which spread all over 
the detector, faked by photomultiplier s noise and 
due to photon scattering. The former, correspon- 
ding to the Signal, produce the Cerenkov photon 
pattern. They are Gaussian distributed, with a 
width a ~ 0.5 cm, reflecting essentially the uncer- 
tainty related to the pixel size. The latter, con- 
stitute an essentially flat background modulated 
by the geometry of the detection plane. Figure 
H shows the distribution of the hit residuals with 
respect to the expected Cerenkov pattern for a 
3 cm thick, 1.03 refractive index aerogel radiator. 
The probability density function for a detected 
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Figure 3. Distribution of the hit residuals with 
respect to the expected Cerenkov pattern. 



hit to belong to the pattern is therefore expressed 
as: 
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Figure 4. Reconstruction of a simulated helium 
event. The reconstructed photon pattern (full 
line) includes both, reflected and non-reflected 
branches. The inner and outer circular lines 
correspond respectively to the upper and lower 
boundaries of the conical mirror. The square is 
the limit of the non-active region. 



P = (1-6) 



exp 



(2) 



where b is the photon background fraction, b/d 
is the background fraction per unit of distance 
(~ 10~ 3 /cm) and n is the closest distance from 
the hit i to the photon pattern. 

Complex photon patterns can occur at the de- 
tector plane due to mirror reflected photons, as 
can be seen on figure [|. The Cerenkov angle re- 
construction procedure relies on the information 
of the particle direction provided by the Tracker, 
which, once extrapolated to the radiator, provides 
an estimation of the mean photon emission ver- 
tex. The tagging of the hits signalling the passage 
of the particle through the solid light guides in the 
detection plane, provides an additional track ele- 
ment. This, can be used as an additional track 
selection criterion. The best value of 9 C will re- 
sult from the maximization of a Likelihood func- 
tion, built as the product of the probabilities that 
the detected hits belong to a given (hypothesis) 



Cerenkov photon pattern ring, 



i h its 



L(9 C 



(3) 



The RICH setup was fully simulated with 
GEANT3 H for radiators with different thick- 
ness and refractive index. In order to trust the 
reconstruction, only events with at least 3 hits 
close (within ~ 1.5cm) to the reconstructed pho- 
ton pattern were selected. Figure [s] shows the 
reconstructed Cerenkov angle for simulated pro- 
tons with various momenta. The Cerenkov an- 
gle single hit resolution, obtained with aerogel 
of 1.03 refractive index and 3 cm thick, was 
A9 C ~ 5 mrad. This resulted in a velocity resolu- 
tion A/3/ (3 slightly better than 0.1% for protons 
(f3 ~ 1) and better than 0.05% for helium nuclei. 
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Figure 5. Cerenkov angle reconstructed as func- 
tion of particle momentum for an aerogel radia- 
tor. 



not to interact in the radiator, p 7 = e~ d ^ z ^l Lint , 
along the radiator thickness and along the pho- 
ton azimuthal angle (<£>). The geometrical accep- 
tance accounts for photons lost through the ra- 
diator walls or totally reflected on media transi- 
tions, photons absorbed by the mirror and pho- 
tons falling into the non-active detection area. It 
is calculated taking into account the portion of 
the visible photon pattern in units of photon azi- 
muthal angle, e geo — Aip/2ir. Figure || shows the 
geometrical acceptance calculated for an aerogel 
radiator of 1.030 and for events within the AMS 
fiducial volume. The extreme variation of e geo 
from event to event is clear. The light guide ef- 
ficiency factor depends on the incidence angle of 
the photons (0 7 ) on the top of the light guide. It 
is calculated using the probability of a given pho- 
ton to get to the photomultiplier cathode once 
it entered the light guide, and by integrating it 
along the reconstructed photon pattern. 



3. Charge reconstruction 

The Cerenkov photons are uniformly emit- 
ted along the particle path inside the dielectric 
medium, L, and their number per unit of energy 
depends on the particle's charge, Z, and velocity, 
/?, and on the refractive index, n, according to the 
expression B: 



dN, 



dE 



7 cxZ 2 L 



1 



1 



(3 2 n 2 



(4) 



Various factors contribute to the loss of some 
of these photons in the RICH; radiator interac- 
tions (e ra d)i geometrical acceptance (s geo ), light 
guide efficiency (ei g ) and photomultiplier quan- 
tum efficiency (e pm t)- Accordingly, the number 
of counted photoelectrons in the detector is given 
by: 



N 7 £ 



rad &geo &lg £pmt 



(5) 



All the efficiency factors, but the PMT effi- 
ciency, depend on the particle direction and of 
its incidence point on the radiator. The radia- 
tor factor depends on the distance, d, traversed 
by the photons inside the radiator. It is calcu- 
lated by integrating the probability of a photon 
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Figure 6. Photon geometrical acceptance for 
events in the RICH detector. 



The charge of the radiating particles is derived 



from the number of photoelectrons, 



close to 



the previously reconstructed photon pattern (see 
section ||). The number of radiated photons is 
obtained by correcting n p . e , by the overall event 
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efficiency, which can be written as: 



4. Conclusions 



Stot 



rad JO 



H,. 



Kpaths 

dz pi 



dcp 



(6) 



where -ffrad is the radiator thickness, n pa ths is 
the number of visible branches constituting the 
reconstructed pattern (i.e. reflected and direct 
branches) and pi being the reflectivity for the i th 
path. 

The charge is then calculated according to ex- 
pression ^, where the normalisation constant can 
be evaluated from a calibrated beam of charged 
particles. In the case of the present results it 
was obtained from 10000 simulated helium nuclei 
events. Figure ^ shows the reconstructed charge 
for elements ranging from Helium to Nitrogen. It 
was obtained using a 3 cm thick, 1.03 refractive 
index aerogel radiator, for events with geomet- 
rical acceptance greater than 60%. The charge 
resolution ranges from AZ/Z ~ 15% for Helium 
to AZ/Z - 5.5% for Nitrogen. 



AMS is a spectrometer designed for anti- 
matter, dark matter searches and for measuring 
relative abundances of nuclei and isotopes. Its 
installation in the International Space Station is 
scheduled to 2005, where it will operate for a 3 
year period. The instrument will be equipped 
with a proximity focusing RICH detector based 
on an aerogel radiator, enabling velocity measure- 
ments with a resolution of about 0.1% and ex- 
tending the charge measurements up to the Iron 
element. The velocity of the cosmic rays is mea- 
sured through the reconstruction of the Cerenkov 
angle using a maximum Likelihood approach. 
The method consists on finding the Cerenkov an- 
gle maximising the overall probability of the de- 
tected hits to belong to its correspon ding pat- 
tern. Charge reconstruction is made in an event- 
by-event basis. It is based both on the velocity 
reconstruction procedure, which provides a recon- 
structed photon pattern, and on a semi-analytical 
calculation of the overall efficiency to detect the 
radiated Cerenkov photons belonging to the re- 
constructed photon ring. 
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Figure 7. Charge reconstruction with simulated 
data in the RICH detector. Gaussian fits are su- 
perimposed to the distributions. 
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